Short-range coherence of a lattice Bose atom gas in the Mott insulating phase 
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We study the short-range coherence of ultracold lattice Bose gases in the Mott insulating phase. 
We calculate the visibility of the interference pattern and the results agree quantitatively with the 
recent experimental measurement (Phys. Rev. Lett. 95, 050404 (2005)). The visibility deviation 
from the inversely linear dependence on the bare on-site interaction Uo is explained both in smaller 
and larger Uo- For a smaller Uo, it comes from a second order correction. For a larger Uo, except 
the breakdown of adiabaticity as analyzed by Gerbier et al, there might be another source to cause 
this deviation, which is the diversity between Uo determined by the single atom Wannier function 
and the effective on site interaction U e s for a multi-occupation per site. 

PACS numbers: 03.75.Lm,67.40.-w,39.25.+k 



The observation of the Mott insulating phase in ul- 
tracold Bose gases in an optical lattice opens a new era 
to investigate exactly controllable strong-correlated sys- 
tems pi Q . For a one-component lattice Bose gases, the 
Bose Hubbard model captures the basic physics of the 
systems 0. The theoretical studies mostly focused on 
the sharp phase transition between the superfluid/Mott- 
insulator [J, IE IE IE IE El ■ This phase transition may 
play an important role in various quantum information 
processing schemes pT|. 

Recently, the residual short-range interference in the 
insulating phase has been predicted by numerical studies 
[l2|. This phase coherence has been observed by a mea- 
surement of the visibility of the interference pattern [l3|. 
It was found that the visibility is inversely proportional 
to the on-site interaction strength Uo of the Bose Hub- 
bard model in a wide range. In explaining their data, 
Gerbier et al assumed a small admixture of particle-hole 
pairs in the ground state of the Mott insulating phase. 
They showed that the visibility of interference pattern 
calculated by this ground state may well match the ex- 
perimental data in a wide intermediate range of Uq. 

There were deviations from the inverse linear power 
law in both small and large U in the measurement of the 
visibility. Gerbier et al interpreted the large U deviation 
is caused by a breakdown of adiabaticity since the ramp- 
ing time used in the experiment has been close to the 
tunnelling time. For the deviation in a small U , there 
was no explanation yet EH • 

In this paper, we will analytically prove the inverse lin- 
ear power law of the visibility for intermediate U in the 
zero temperature. Here the words 'intermediate V (as 
well as 'small U\ 'large IP in this work) mean the mag- 
nitude of U — U c is intermediate (small or large) , with U c 
the critical interaction strength of the superfluid/Mott 
insulator transition. The result is exactly the same as 
that obtained by Gerbier et al by assuming a small ad- 
mixture of the particle hole pair in the ground state ^| . 
We also show the deviation of the visibility from the in- 
verse linear power law in a small U is caused by a second 
order correction. For the large U, we show that, except 
the explanation by the authors of the experimental work, 



owing to the multi-occupation per site, the effective on- 
site interaction Ueg which appears in the Bose Hubbard 
model \13L lltx Il7| is different from Uq which was deter- 
mined by the single atom Wannier function and used to 
fit the data of the experiment. 

We consider a one-component Bose gas in a 3- 
dimensional optical lattice described by a periodic po- 
tential Vo(r). Although the real experimental system 
was confined by a trap potential, we here only pay our 
attention to the homogeneous system. Beginning with 
the expansion of the boson field operators in a set of 
localized basis, i.e., tp(r) — J2i a i w (^ — ^i) an d keep- 
ing only the lowest vibrational state, one can define 
an on-site free energy / = hi + Uh(h — l)/2, where 
n is the average occupation per site. The on-site en- 
ergy / and the bare on-site interaction U are defined by 
I = j drw*{f)[-^V 2 + V (f)]v)(?), 



r = tEBlL I dr >(f)| 4 



(i) 



This on-site free energy contributes to the chemical po- 
tential by fi = —df/dh and defines the effective on-site 
interaction HE EE El 



U, 



eff 



d 2 f/dn 2 



(2) 



For the single occupation per site, U c s = U = Uq and 
the difference appears for n > 1. We will be back to this 
issue later. The Bose Hubbard model for a homogeneous 
lattice gases is defined by the following Hamiltonian 



H 



+ \ " t i Uctf \ - + + 



(3) 



where (ij) denotes the sum over the nearest neighbor 
sites and \x is the chemical potential. The tunnelling 
amplitude is defined by 



tB,F; 



dfw*(r + fi)[- 



h 2 v 2 

2m 



V (r)]w(r + r.j), 



for a pair of the nearest neighbor sites 



2 



Our main goal is to calculate the interference pattern 



(4) 



which is related to the density distribution of the expand- 
ing atom clouds by p(r) — j^—\w(k = mr/ht ex )\ 2 S{k) 
with to the atom mass and t ex the time of the atom free 
expansion [T2I Il8l | . Since we are interested in the Mott 
insulating phase, we can calculate S(k) by taking the tun- 
nelling term as a perturbation. To do this, we introduce 
a Hubbard-Stratonovich field in the partition function 

Z[J,J*} = J V<S>*V<f>Va*T>aexp^-S 

Jo do) Jo 

-tj o dT(a;-$t + Ji/t)fa-*j + Jj/t)\, (5) 

where So is the t-independent part in the full action and 
J and J* are currents introduced to calculate correlation 
functions. Integrating away at and a* and transferring 
into the lattice wave vector and thermal frequency space, 
one has 

Z[J*,J] = |p$*P$exp|^(-$i n e- 1 (fc, i c„)$,- n 



4- T* + /- A* -I 72 7- 

+ J fc,n fc,n + J k,n V k,n + £fc J k,n J k,n> 



(6) 



where tk = — 2t ^2 V=X y z cos The correlation function 
is calculated in a standard way: 



(a~ ar ) — 

s k,n k,nl 



1 (5 2 Z[J*,J] 



Z[0,0] 
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-0(*,*w„) + — . (7) 



The interference pattern then may be expressed as 



(8) 



In the Mott insulating phase, the correlation function 
Q(k, iuu n ) has been calculated by slave particle techniques 



(9) 



Q 1 (k,jw n ) = e k -£kV(a + l)T 
^ — ' 7 



a=0 



iliJn + jJL — all 



where the slave particle occupation number is given by 

1 



exp{/3[-iA - a/i + a(a - l)U c s/2]} ± 1 ' 



(10) 



which obeys n a = 1 and ^ a an Q = N in the mean 
field approximation |19| . A is a Lagrangian multiplier 



to ensure J2 a ™ Q = 1 ■ The sign ± corresponds to the 
slave fermion or boson, respectively. In previous works, 
we have show that the slave fermion approach may have 
some advantages to the slave boson approach^ [l(|. We 
then take the slave fermion formalism. In the Mott insu- 
lating phase, since U e g, [i ;§> t, one can expand G{k, iu n ) 
in terms of ek/(ito n — (J- + 0-U e s) and the interference pat- 
tern reads 



S(k) 



-I u,„)-i) 
1^^(-1)»^(AK))» +1 , (11) 



A{u> n ) 



a=0 



(<*+!)- 



/i - aU cS ' 



Making the frequency sum, one has, to the first order 
of efc, 



S(k) « - ^ n B (aU e s -(*)((* + l)(r 

a 

- e fc /3^[( a + l) 2 K+ 1 -n Q ) 1 



n a ) 



x riB(aU e H - + riB(aUcS - ^)) (12) 
2e fc \ - 



TTZ z2( n B(aU cff - (J,)- n B {~lU cS - /Lt)) 



a<7 

x (a+l)( 7 + l)(n a+1 -n a )(n 7+1 -n 7 ), 

where n B (aU cS - //) = [e^ ^"^) - l]" 1 . In the limit 
T — > and the n -th Mott lobe, one knows (n — l)U e g < 
fi < noU c s and n a — <5 a ,„ . Substituting these into H12II . 
one obtains the zero temperature value of S(k) 

S(k 1 T = 0) = n -2n (n + l)^ . (13) 

U eff 

This is what Gerbier et al obtained by assuming the 
particle- hole pair admixture in the ground state |13| . In- 
tegrating along one lattice direction, the corresponding 
2D visibility is given by 



V 



Pn 



S n 



S n 



lf4 (14) 

U off 



Pmax ~r Pmin &m 

for z = 6, where p max and p mm are chosen such that 
the Wannier envelop was cancelled. This is the inverse 
linear power law used to fit the experimental data |13| . 
However, the experimental data deviated from this power 
law fit when U e s/zt < 8. In terms of (|1 1|) . we think 
that this comes from a second order correction. A direct 
calculation shows that the second order correction in zero 
temperature is given by [2£j 

S^S(k) =3no(n + l) 2 J-. (15) 

Thus, the 2D visibility for no = 1 is modified to 



3[/- eff (l + 32[/- s <73)' 
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with U c ff = Ueg/zt. In Fig. ^ we show the visibility 
against U c g in a log- log plot for n = 1. This second 
order correction suppresses the visibility for a small U c s 
while the exponent of the power law seems deviating from 
— la little. These features agree with the experimentally 
measured data. 



0.4 



> 




FIG. 1: Visibility of the interference pattern versus U e g ac- 
cording to 11611 in a log-log plot (the dot line with circles) . 
The solid line is the inverse linear power law ()14fl and the 
dash line is a power law fit with an exponent —0.95 to 11611 . 

We have neglected the finite temperature effect to com- 
pare with the experiment although our theory is in finite 
temperature. In fact, there may be a finite temperature 
correction to the interference pattern in the second order. 
According to (fTO . it is given by, near n = 1 

6^S T (k) = lS^V (17) 

which may further suppresses the visibility. For instance, 
at T = l.Ozt - lf^nK, the ratio between ([TT| and lO 
is 

S T )S(k ) = 3(n x ) V/2 

8<!*)S(k) 

= 0.106,0.098,0.087, and 0.064 

for C/ofi = 6, 7, 8, and 10. However, the temperature in 
the Mott insulator is difficult to be estimated in the ex- 
periment |22|. Thus, a quantitative comparison of the 
finite temperature calculation to the experiment data is 
waiting for more experimental developments. 

We now discuss the large U deviation from the inverse 
linear power law, which has been seen in the experiment 
and explained by the breakdown of adiabaticity ^1 • We 
will reveal another possible source for this deviation. As 
we have mentioned before, the value of £/ e ff may be dif- 
ferent from U and Uq for n > 1. Our above calculation 
showed an inverse linear power law to U e g whereas the 
experimentalists used Uq to fit their data. 




AVERAGE OCCUPATION 

FIG. 2: (Color online) The effective on-site interaction U e a 
versus the average occupation per site, n in a n-log(f/ e fr) 
plot. The thin solid lines are linear fits to variational data for 
V = 11.95, 14.32, 16.25 and 29 E R (empty circles, filled trian- 
gles, empty triangles and filled circles, respectively). The dash 
line is critical interaction strength calculated by the mean field 
theory |j, |5J . The thick horizontal lines are the on-site inter- 
actions Uq calculated by the single atom Warrier function. 



Due to the interaction, the atom energy band may be 
modified and the Wannier function may be broadened, 
compared to the single atom ones. In Ref. 01 > we 
have considered the mean field interaction and made a 
variational calculation to the Wannier function by using 
Kohn's method [2l| ■ The direct result of the broadening 
of the Wannier function is the bare on-site interaction 
U becomes weaker than Uq which is calculated by the 
single atom Wannier function. The n-dependence of I 
may further reduce U e g from U. In Fig. [21 we plot 
U c g versus n. In the low part of Fig. [2J three typi- 
cal lattice depths are considered, Vb = 11.95,14.32 and 
16.25 Er{= %^), corresponding to the critical interac- 
tion strengths of the uq = 1,2 and 3 Mott states. The 
up-part is for Vq = 29Er, which was the lattice depth 
where the adiabaticity breaks 0] . 

Several points may be seen from Fig. [21 First, the crit- 
ical values of Vq — IA.32Er for n — 2 and 16.25Er for 
uq — 3 are closer to experimental ones, 14.1(8) En and 
16.6(9) E R [13, comparing to 14.7 E R and 15.9-Er, corre- 
sponding to the single atom Wannier functions. Second, 
the variational data are downward as n indicates that 
— logC/cff > — logC/o for n > 1. This may cause two re- 
sults: (a) If — log U c ff deviates from — log Uq a small mag- 
nitude, the power law fit presents an exponents —(1 — 5). 
This has been observed in experiment, which is —0.98(7) 
[l3j . (b) As n increases, the deviation becomes signifi- 
cant. This may appear in a large Vq. In the experiment, 
the latter appeared in Vq > 29Er. We show that, in Fig. 
[3 the deviation is not a small magnitude for Vq = 2%Er. 

In summary, we studied the short-range coherence in 
the Mott insulating phase with a finite on-site interaction 
strength. The interference pattern and then its visibility 
were calculated by using a perturbation theory. The in- 
verse linear power law of the visibility to the interaction 
strength, which was found in the experiment, was exactly 
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recovered. We further discussed the deviation from this 
power law both in a small and large Uo. We found that 
a second order effect suppresses the visibility for a small 
Uq while its up-deviation in a large Uq might be caused 
by the difference between Uq and U e s except the possible 



breakdown of adiabaticity. 
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